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Abstract. In this letter we present an idea which reconciles a 
homogeneous and isotropic Friedmann universe with a fractal 
distribution of galaxies. We use two observational facts: The 
flat rotation curves of galaxies and the (still debated) fractal 
distribution of galaxies with fractal dimension D = 2. Our idea 
can also be interpreted as a redefinition of the notion of bias. 
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It is known since twenty years that the galaxy distribu- 
tion exhibits fractal behavior on small scales (Peebles 1980, 
Mandelbrot, 1982). Several recent statistical analyses of three 
dimensional galaxy catalogues indicate that galaxies are dis- 
tributed fractally with dimension D ~ 2 out to the largest 
scales for which statistically significant data is available, i.e., 
up to about 100-200/i _1 Mpc (Coleman & Pietronero, 1992; 
Sylos Labini et ai, 1998). In the opinion of the authors of the 
present note, no data is pointing convincingly towards a ho- 
mogenization of the galaxy distribution ( i.e., D — 3 ). We 
believe that the standard analyses which indicate that fluctu- 
ations of the galaxy number density decrease on large enough 
scales are not well suited to study scale invariant structures, 
since they a priori assume the existence of a well defined av- 
erage density inside the given sample (Peebles, 1993; Davis 
1997). On this point the scientific community has not reached 
consensus (Wu et ai, 1998). In this work we assume without 
further arguments that the fractal picture is correct. Our main 
point here is to show that a fractal distribution of galaxies even 
up to the Hubble scale may be consistent with a homogeneous 
and isotropic universe. 

On the other hand, many observations, most notably the 
superb isotropy of the cosmic microwave background together 
with its perfect blackbody spectrum, give strong evidence 
that the geometry of our Universe is very homogeneous and 
isotropic on large enough scales, a so-called Friedmann model. 

In this letter we argue that this seemingly flagrant contra- 
diction of different pieces of observational data may actually 
be reconciled in a rather simple way. 

Our main point is that a fractal distribution of galaxies 
need not imply a fractal matter distribution (as it is also 
pointed out in Wu et at, 1998). 

It is often mentioned that galaxies represent peaks in the 
matter distribution. Let us compare these to the distribution of 
mountain peaks on the surface of the earth which we know are 



fractally distributed over a certain range of scales. But it would 
be false to conclude from this fact that the radius of the surface 
of the earth from its center is grossly variable; as we know, this 
number is very well approximated by a constant. A similar 
effect may actually be at work in the matter distribution of 
the universe. 

We first interpret the COBE DMR results (Bennett et 
ai, 1996) and other observations of CMB anisotropies on 
smaller angular scales (De Bernardis et al, 1997). They indi- 
cate that matter fluctuations are reasonably well described by 
a Harrison- Zel'dovich spectrum with (5M/M)(X s (t)) ~ 10" 4 , 
where \n(t) denotes the comoving horizon scale, Ah(£) = r\ ~ 
t/a, and a is the scale factor. But small density perturba- 
tions in a Friedmann universe grow only once the universe 
becomes matter dominated and even then rather slowly (pro- 
portional to the scale factor a tx t 2 ^ 3 ). Density fluctuations 
today on a given comoving scale A should thus be on the or- 
der of \- [i/2 6(\) = 10- 4 (toA(A)) 2/3 ~ 10- 4 (t7o/A) 2 . Here t(X) 
denotes the cosmic time at which the scale A enters the hori- 
zon; to and r/o are the present cosmic and conformal times 
respectively. For example on the scale A = 200/i~ Mpc matter 
density fluctuations should not be larger than 
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A=200h- 1 Mpc 



^ U=200h- 1 Mpc 

(3000h~ 1 Mpc/200/T 1 Mpc) 2 10~ 
0.015 



today. (Here h parameterizes the uncertainty of relating the re- 
cession velocity of a galaxy to its distance, i.e. the uncertainty 
in the Hubble parameter Ha = 100/ikm/s/Mpc.) 

We now argue that a fractal galaxy distribution may well 
be compatible with such a homogeneous Universe. This can be 
seen by the following observation: It is well known that the 
dark matter distribution around a galaxy leads to flat rotation 
curves (Rubin et al., 1980). Without cutoffs this yields a matter 
density distribution 

Phaio(r) oc r~ 2 

around each galaxy. To obtain the total matter distribu- 
tion, we actually have to convolve phaio with the number dis- 
tribution of galaxies, tig- The fractal dimension D = 2 means 
that the average number of galaxies in a sphere of radius r 
around a given galaxy, denoted by Nq(t), scales like r 2 . This 
implies that the mean number density of galaxies around an 
occupied point decays like 
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Our main finding is the simple fact that the convolution of 
these two densities gives a constant (up to logarithmic correc- 
tions) 



p = 71G * Phaio oc 1 



(1) 



More precisely, for no,(r) = C/r and phalo(r) = A/r 2 we 
obtain ( with |y| = y and |x| = a;) 



p(x) =AC id 
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= 4nAC[l + ln(7? max /x)] (2) 



This shows that a fractal galaxy distribution with dimension 
D — 2 together with flat rotation curves, indicate a smooth 
matter distribution in agreement with our expectations from a 
Friedmann universe with small fluctuations. The dark matter 
distribution of a two dimensional model where the galaxies are 
distributed with fractal dimension D = 1 is shown in Fig. |l| 
Clearly, this dark matter distribution is very homogeneous up 
to finite size effects. 

Note that the essential ingredient for this result is that 
riG oc l/r a and Phaio(?") oc 1/r 13 with a + (3 — 3. 



a sufficient distance from a given galaxy the only measurable 
density is the total density which contains relevant contribu- 
tions from many galaxies. But in principle it is possible to 
assign to each galaxy a density profile phaio, and it is interest- 
ing to note that the form of the density profile indicated by 
measurements, which are possible close to isolated galaxies, is 
just such as to lead to a constant total density if we convolve 
it with the number density of galaxies. 

Another objection may be, that with this density profile 
the total mass of a galaxy is infinite. But this is in fact ir- 
relevant, since not only the mass of a single galaxy goes to 
infinity as r — » oo, but also the galaxy density goes to zero 
as r — > oo and this just in a way that the measurable total 
density is constant. Besides, the integral of phaio should not be 
considered as the 'mass of the galaxy'. More profoundly, the 
flat rotation curves are a consequence of the fact that the grav- 
itational potential remains constant during linear clustering in 
a Friedmann universe, which holds beyond the scale of single 
galaxies. 

Clearly, the amplitudes C and A of na and phaio depend 
on the type of galaxies considered. Galaxies of different abso- 
lute luminosities in general have different circular speeds and 
different abundances. 

To quantify Eq. (^) we use the Tully-Fisher relation (Tully 
& Fisher, 1977), between the circular speed v of a galaxy and 
its luminosity L, v = v»(L/L*)°- 2B , with v, = 220km/s. The 
luminosity L» corresponds to an absolute magnitude of M* ~ 
-19.5. 

The abundance of L» galaxies as estimated from Sylos 
Labini et aZ. (1998) is 



n 

Na(r),=By with £>* ~ 2 



Considering for the time being just L« galaxies, this gives 
n G = C,/r with C» ~ 3B„/47r. 

To determine the amplitude of the halo density, phaio = 
A/r 2 , we use the flatness of galaxy rotation curves with (Rubin 
et al, 1976) 



Uv t /c) 2 ~ 10" 7 = GM{r)/r 



Fig. 1. We show the dark matter distribution of a two dimensional 
set of fractally distributed galaxies (D = 1,uq oc 1/r) (filled cir- 
cles) each surrounded by a dark matter (dots) halo with distribution 
Phaio 1/r. The halos sum up to a very homogeneous dark matter 
distribution. 

It is important to insist that the cutoff i? max is larger than 
all the scales x considered. Far away from galaxies, for exam- 
ple in a void, many galaxies contribute to the density in a 
given point and the important message is just that flat rota- 
tion curves indicate actually that the total resulting density 
may be rather constant in voids and has about the same value 
as it has close to galaxies. 

One may argue that the notion of a certain lump of matter 
'belonging' to a certain galaxy is ill defined sufficiently far away 
from a galaxy and thus the 'halo' density cannot be defined on 
scales larger than, say, half the distance to the next galaxy. 
With this objection we agree in practice, it just means that at 



where c denotes the speed of light. Combining this with 
M (r) — 4nAr gives 



A t = 



K/c) 2 
8ttG 



The convolution of hg with phaio then leads to 
3B„(^/c) 2 



p = n G * p halo ~ 47rC*A* 



8ttG 



(3) 



Comparing this with the critical density of a universe ex- 
panding with Hubble parameter Ho, 



Pc = 



3# 2 c 2 



(4) 



8nG ' 

we obtain a density parameter of order unity, 
n = p/p c ~ vlB,/Hl ~ 1 . (5) 
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Clearly, this estimate is very crude since not all galaxies 
have the same rotation speeds and the constant B depends on 
the luminosity of the galaxy. But it is a reassuring non-trivial 
'coincidence' that the density parameter obtained in this way is 
of the right order of magnitude. To refine this model we have 
to find a Tully-Fisher type relation for B(L) and integrate 
over luminosities. But since the abundance of galaxies decays 
exponentially with luminosity above L*, we have 



v(L) 2 dB(L) ~ vlB, 



and reproduce the result (g). 

Our arguments suggest that a fractal galaxy distribution 
may well be in agreement with a smooth matter distribution. 
Neglecting log corrections (which may be absent in a more re- 
alistic, detailed model and which are certainly not measurable 
with present accuracy), our model describes a perfectly homo- 
geneous and isotropic universe. We do not specify the process 
which has induced small initial fluctuations and finally led to 
the formation of galaxies. We are thus still lacking a specific 
picture of how the fractal distribution of galaxies may have 
emerged. Purely Gaussian initial fluctuations are probably not 
suited to reproduce a fractal galaxy distribution. But cosmic 
strings or other 'seeds' with long range correlation could do it. 
A working model remains to be worked out. 

It may be useful to mention that the view presented here 
actually redefines the notion of bias. In the standard scenario, 
the dark matter density field (on large scales) is Gaussian and 
galaxies form in the peaks of the underlying distribution and 
their correlation functions are simply relatedlO Here we con- 
sider the possibility that the galaxy and dark matter distribu- 
tions may have different correlation properties and hence dif- 
ferent fractal dimensions, namely D — 2 for the galaxies and 
D = 3 for dark matter. Especially in view of Eq. (hi) , we want 
to warn the reader against interpreting the galaxy distribution 
as proportional to the matter distribution even on large scales. 

A similar idea is the one of a universe with a fractal galaxy 
distribution but a dominant cosmological constant A > 8irGp. 
This possibility is also discussed in Baryshev et al. (1998), in 
connection with the linearity of the Hubble law. 

More precisely, the fractal dimension of a set of density 
fluctuations can depend on the threshold (see Fig. |^). In a 
realistic model, we would expect that also the dark matter, 
above a certain threshold is fractally distributed. In galaxy 
catalogues, this tendency is actually indicated. Observations 
show a slight increase of the fractal dimension with decreasing 
absolute luminosity of the galaxies in the sample (Sylos Labini 
et al., 1998), however, still with relatively modest statistics. 
Such a scenario is naturally formulated within the framework 
of multi- fractals (Falconer, 1990; Sylos Labini et al., 1998). 
It is, for example, a well known fact that bright ellipticals lie 
preferentially in clusters whereas spiral galaxies prefer the field 
(Dressier, 1984). From the perspective of multi-fractality this 
implies that ellipticals are more clustered than spirals, i.e., 
their fractal dimension is lower (Giovanelli et al., 1986). 

Our arguments indicate that the voids might be filled by 
dark matter. But since this dark matter is relatively smoothly 
distributed, it cannot be detected by measurements sensible 
only to density gradients (like, e.g., peculiar velocities). One 
needs to determine the total density of the universe, for exam- 
ple by measuring the deceleration parameter go- 



10 



£10 



10' 



D=1.1 
D=1.4 
D=1.6 
D=1.8 
D=1.9 
D=2 



10 - 



10" 



10" 



Fig. 2. The mean density around an occupied point (T(r)) is shown 
for our two dimensional model. The galaxy density with fractal di- 
mension D = 1 is shown as filled circles. As the density threshold 
decreases the fractal dimension approaches D = 2. The chosen over- 
densities for the fractal demensions of 1.1, 1.4, 1.6, 1.8, 1.9 and 2 
are 2.1, 1.76, 1.55, 1.34, 1.22 and 1 respectively. 
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